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Abstract 


The Mars Science Laboratory (MSL) Entry, Descent, and Landing Instrumentation (MEDLI)/Mars 
Entry Atmospheric Data System (MEADS) project installed seven pressure ports through the MSL 
Phenolic Impregnated Carbon Ablator (PICA) heatshield to measure heatshield surface pressures 
during entry. These measured surface pressures are used to generate estimates of atmospheric quan- 
tities based on modeled surface pressure distributions. In particular, the quantities to be estimated 
from the MEADS pressure measurements include the dynamic pressure, angle of attack, and angle 
of sideslip. This report describes the calibration of the pressure transducers utilized to reconstruct 
the atmospheric data and associated uncertainty models, pressure modeling and uncertainty anal- 
ysis, and system performance results. The results indicate that the MEADS pressure measurement 
system hardware meets the project requirements. 
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Chapter 1 

Introduction 


This report documents transducer calibration, system error modeling, pressure, pressure uncer- 
tainty models and system performance results and uncertainty quantification for the Mars Science 
Laboratory (MSL) Mars Entry Atmospheric Data System (MEADS) [1,2]. MEADS is a seven- 
port flush air data system (FADS) integrated with the 4.5 m diameter MSL aeroshell, the largest 
aeroshell ever flown to Mars [3]. The Mars Science Laboratory entry vehicle was comprised of a 
70 degree sphere-cone heatshield and backshell consisting of a stack of three truncated cones. The 
forebody was similar to the heatshield geometry developed for the Viking Mars landers. Phenolic 
Impregnated Carbon Ablator (PICA) was used for the Thermal Protection System (TPS) mate- 
rial. The backshell configuration was also similar to Viking, with a third cone section added to 
accommodate the parachute volume. The MSL vehicle as-built outer mold line is shown in Fig. 1.1. 
During most of entry, the capsule used a radial center of mass offset to fly at an angle of attack 
(approximately 16 degrees at hypersonic conditions). This attitude produced lift to fly a guided 
entry profile, reducing the landing footprint to a much smaller size than any previous Mars mission. 
To fly the guided entry, the vehicle carried four pairs of reaction control system jets to perform 
maneuvers and damp rates. The four pairs of jets could be fired rapidly in different combinations 
to provide control torques about roll, pitch, yaw, or any other axis by modulating the pulses of the 
jet. 

The MSL heatshield differs from previous Mars mission designs in that the TPS is based upon 
a tiled scheme, bonded directly to the aeroshell structure with a room temperature vulcanization 
(RTV) sealant. Viewed exteriorly, MEADS surface pressure orifices are holes of approximately 2 
mm diameter drilled directly through the PICA TPS tiles. Metallic spools through the carbon 
fiber-on-aluminum honeycomb aeroshell structure complete the flowpath from surface orifice to 
transducer tube as shown in Figure 1.2(a). MEADS therefore represents somewhat of a departure 
from traditional aircraft FADS, or even the Shuttle Entry Air Data System (SEADS) [4-7], where 
orifices are precisely machined to relatively small diameters in smooth metallic or nonablating 
materials with negligible surface roughness. The MEADS ports are arranged in a cruciform pattern 
over the heatshield windside (see Figure 1.2(b)), in locations where the boundary layer is predicted 
to remain laminar and aeroheating is predicted to remain relatively low over the expected entry 
trajectory [1]. All of the pressure ports are located a minimum of 3-inches from the PICA tile 
seams. 

Figure 1.3 shows a timeline of the different events along MSL’s entry, descent, and landing 
(EDL) trajectory [8]. Note that the event times correspond to the nominal trajectory. Actual 
times differed from the nominal event times due to dispersions in the entry time and atmospheric 
uncertainties. EDL consists of six major segments: Exo- Atmospheric, Entry, Parachute Descent, 
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Spacecraft Coordinate Frame 
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(a) MEDLI/MEADS Tube Configuration (b) MEDLI/MEADS Pressure Port Locations 

Figure 1.2: MEDLI/MEADS Geometry 


Powered Descent, Sky Crane, and Fly Away [9]. The Exo- Atmospheric segment begins once the 
cruise stage separation command is sent. Once the cruise stage separates, Guidance, Navigation, 
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and Control (GNC) is enabled. Once enabled, the entry body is despun and turned to its entry 
attitude. Then, the two 75-kg Cruise Balance Masses (CBMs) are jettisoned to enable aerodynamic 
lift. 
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(a) Cruise Stage Separation to Parachute Deployment (b) Parachute Deployment to Touchdown 

Figure 1.3: MSL Entry, Descent, and Landing 


The Entry segment starts with the vehicle at the Entry Interface Point (EIP) defined at 3522.2 
km from the center of Mars, approximately 540 seconds after cruise stage separation. During the 
Entry segment, the vehicle goes through peak heating and peak deceleration, the Reaction Control 
System (RCS) controls the lift vector to achieve the desired down-range and cross-range target. 
Just prior to parachute deployment, six 25-kg Entry Balance Masses (EBMs) are jettisoned to 
eliminate lift and the vehicle rolls to point the Terminal Descent Sensor (TDS) to the ground. This 
maneuver is called the Straighten Up and Fly Right (SUFR) maneuver. Figures 1.4 and 1.5 shows 
the nominal trajectory over the time period from entry interface to parachute deployment. 

The Parachute descent segment starts with the parachute deployment triggered once the vehicle 
reached Mach 1.7. Once the vehicle achieves a speed of Mach 0.7, the heat shield is jettisoned and 
the TDS starts acquiring the ground. Note that the MEDLI instruments are powered off 10 seconds 
prior to heat shield jettison. The command to jettison the backshell and the parachute is issued at 
an altitude of 1.6 km and at a velocity of approximately 79 m/s. Just before backshell separation, 
the Mars Landing Engines (MLEs) are primed in preparation for the start of the powered descent 
segment. 

The Powered Descent segment begins at backshell separation. During powered descent, eight 
independently throttleable MLEs are actuated, initially to execute a divert maneuver for backshell 
avoidance which brings the vehicle to vertical flight at a descent rate of 32 m/s. Once vertical 
flight is achieved, a descent at constant velocity to adjust for altitude error at backshell separation 
starts. This constant velocity phase is followed by a constant deceleration phase, which reducs the 
vehicle’s speed to 0.75 m/s in preparation for the sky crane segment. At this time, the four inboard 
MLEs are throttled down to near shutdown (1%) while the four remaining MLEs were throttled 
at 50%. The Sky Crane segment starts following issuance of the rover separation command, which 
occurrs at an altitude of approximately 18.6 m. The rover is lowered to 7.5 m below the descent 
stage. Then, the descent stage continues to descend until post-touchdown is detected. 

The Fly Away segment starts after touchdown is sensed. Once the descent stage stops it’s 
vertical motion, the bridle and electrical umbilical devices are cut and two of the MLE engines are 
throttled up to 100% while the other two engines are at slightly less than 100%. This causes the 
descent stage to pitch to 45 degrees. Once the turn maneuver is completed, all four engines are 
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throttled up to 100%. Constant thrust is applied to ensure the descent stage impacts the surface 
at least 300 nr from the landing point. 

MEADS collects pressure data at a rate of 8 Hz over the portion of MSLs entry from atmo- 
spheric interface to 10 seconds before heatshield separation, the latter which occurs at a freestream 
Mach number of approximately 1.5. As a basic objective, the MEADS hardware is intended to 
be capable of furnishing data that, without any additional data sources, can be used to recover 
aerodynamic angles and dynamic pressure. Accuracy to within 0.5 degrees in angle of attack and 
sideslip are required over portions of the trajectory where freestream dynamic pressures exceed 
1750 Pa and 1250 Pa, respectively. Dynamic pressure measurements are required to within 2% 
where the freestream value exceeds 850 Pa. Used in combination with inertial measurement unit 
(IMU) data and analyzed with Extended Kalman Filter (EKF) methods [10-14], MEADS data will 
be used for post-flight reconstruction of vehicle attitude, dynamic pressure, Mach number, atmo- 
spheric density, and wind velocity [15]. As MSL is breaking new ground as the first tiled ablator 
NASA has flown on an entry vehicle heatshield, MEADS is unique in its distinction as the first 
heatshield-integrated FADS to be employed on an extraplanetary mission on an aeroshell larger 
than and in an aerothermal environment more extreme than experienced in any previous Mars 
mission [16]. 

This document describes how the uncertainties in the as-built MEADS hardware system and 
pressure models contribute to the uncertainties in the reconstructed data. The remainder of this 
document is organized as follows. Chapter 2 describes the basic MEADS data processing algorithm, 
in which pressure measurements are used in conjunction with a model of the surface pressure 
distribution (based on CFD models) to determine the freestream atmospheric conditions (dynamic 
pressure, static pressure, angle of attack, and angle of sideslip) using a nonlinear least-squares 
algorithm. Also in Chapter 2, the CFD pressure models are described along with some description 
of experimental validation. 

The MEADS pressure transducer calibration is described in Chapter 3. A basic overview of 
the calibration methodology is described, and results from thermal vacuum testing are given along 
with an uncertainty analysis. Chapter 4 describes hardware system error modeling and CFD 
pressure distribution uncertainty models. These uncertainty models are then utilized in Chapter 
4 to quantify how these uncertainties map to uncertainties in the reconstructed air data state. 
Concluding remarks are given in Chapter 5. 
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Figure 1.4: MSL Nominal Trajectory: Altitude and Mach Number 
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Figure 1.5: MSL Nominal Trajectory: Dynamic Pressure and Aerodynamic Angles 
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Chapter 2 

Atmospheric State Estimation and 
Surface Pressure Modeling 


The basic MEADS atmospheric state estimation procedure is essentially similar to the Shuttle 
Entry Air Data System (SEADS) state estimation described in [4], a point-wise least squares fit of 
the surface pressure measurements to the modeled pressure distribution to determine minimum- 
variance estimates of the atmospheric states. This algorithm requires a model of the surface pressure 
distribution at each port location as a function of the free stream flight conditions. For MSL, this 
model consists of tabulated CFD pressure solutions for a range of flight conditions that encapsulate 
the expected trajectories. An uncertainty model of the pressure distribution and sensor hardware 
is also required for appropriate weighting of the measurements in the least-squares algorithm. 

The core of the MEADS pressure model is based on CFD predictions run at representative 
nominal Mars entry conditions, while the pressure uncertainty model borrows from off-nominal 
Mars entry CFD, wind tunnel experimentation, and CFD simulation of wind tunnel experiments. 
Surface pressures are extracted from CFD solutions at the MEADS port locations, which are well 
known from flight article laser scan metrology. Comparisons between experimental and numerical 
results were to be used to validate computations and model estimates of error in CFD simulations 
at flight conditions, however these results are beyond the scope of this report. 


2.1 Least— Squares State Variable Estimation 

By defining the atmospheric state vector as x = [a, (3, qoo,Poo] T , P as the vector of N observed 
surface pressures, and the pressure model as h(x,t) then the pressure measurement model can be 
written compactly as 

p = h(x,t k ) + e (2.1) 

where e is the vector of pressure measurement errors and t k is the measurement sampling time 
at time increment k. Definition of the aerodynamic angles is shown in Figure 2.1. The pressure 
measurement model can be approximated by means of the truncated series expansion 

p « h(x,t k ) + H(t k ) (x~x) + e (2.2) 

where x is a reference state and H(t k ) is the Jacobian matrix 

(2.3) 

X=x 


H(t k ) = 


dh 

dx 
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ac 


The problem can be reduced to that of a linear regression problem 

y = Hx + e (2.4) 

where y = p — h (x,tk) + Hx. By virtue of the Gauss-Markov theorem, the best linear unbiased 
estimate of x is the weighted least-squares solution [17], 

x=[H T R 1 Hy l ^R^y (2.5) 

where R is the pressure measurement error covariance matrix. Since the original relationship 
between the states and the measurement in Eq. (2.1) is nonlinear, the estimation scheme can be 
iterated until convergence by successively replacing x by x, using the converged state estimate 
from the previous sampling time (ifc-i) as the initial reference state x. The state estimate error 
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covariance matrix P can then be computed from 

P=(H T BT 1 Hy l ( 2 . 6 ) 

Note that the atmospheric state vector formulation is not unique. In particular, any two vari- 
ables from the set of poo-, pt, q oo, 77, and can be used in the estimator in addition to a and (5. 
The remaining state variables not used directly in the estimator can be calculated as parameters, 
having first estimated the atmospheric state variables. Note that the SEADS algorithm [4] made 
use of Poo and pt in the estimator formulation. 


2.2 Nominal Flight CFD Pressure Modeling 

The MEADS pressure model is built on CFD solutions from the MSL nominal aerodatabase, com- 
puted at points of interest along reference trajectory 05-22 as shown in Fig. 2.2(a). Computations 
were carried out with the Langley Upwind Relaxation Algorithm (LAURA) code [18] using a two- 
temperature thermo chemical nonequilibrium gas model representative of the Martian atmosphere 
with a freestream composed of 97 percent CO 2 and 3 percent N 2 by mass. A representative 
LAURA solution is shown in Fig. 2.2(b). A total of 8 species (CO 2 , CO, N 2 , O 2 , NO, C, N, O) 
were allowable to model the chemical reactions which take place at energetic flow conditions [19]. 
A super-catalytic wall boundary condition was used, which imposes recombination of species to 
freestream composition at the wall so as to yield conservative aeroheating predictions [20]. The vi- 
brational temperature relaxation model of Carnac [19] was used throughout, as it has been assumed 
more accurate than the Millikan- White model [19], however analysis validating differences between 
the models is beyond the scope of this work. Thin-layer Navier-Stokes (TLNS) equations were used 
in the computations, a subset of the full Navier-Stokes equations which omits viscous cross terms 
and thereby reduces the number of calculations required to converge a solution. Previous work had 
shown that differences between solutions computed with the TLNS and full Navier-Stokes equations 
were insignificant for the purpose of computing aerodynamics [19]. During a LAURA computa- 
tion, grids are shock-fitted to the evolving solution with an iterative grid adaptation scheme which 
maintains orthogonality and enforces Re ce u ~ 0(1) at the wall for aeroheating computations. Grid 
adaptation cycles are repeated until integrated aerodynamic coefficients are constant to within 
three significant figures, after which final numerical convergence is obtained using a line relaxation 
algorithm [21], The nominal aerodatabase runs were computed for the Outer Mold Line (OML) 
on a 7-block volume grid with 0.23 million points. During design iterations leading up to the final 
flight configuration (designated OML13-F), the flight TPS was changed from a monolithic ablator 
design to tiled Phenolic Impregnated Carbon Ablator (PICA) design, changing OML proportions 
slightly. Point check cases were run at nominal 05-22 trajectory conditions on a volume grid more 
closely describing the flight article heatshield (MSL-OML13F-21806-clean-heatshield.g) with a to- 
tal of 1.25 million points. Fig. 2.3(a) and (b) show comparisons between point check and nominal 
aerodatabase CFD results are used to obtain an estimate of MEADS pressure sensitivity to changes 
in the flight article OML, grid topology, and grid refinement. 

Numerical solutions at off-nominal flight conditions were used to bound the expected ranges of 
MEADS port environments and sensitivities. LAURA solutions from runs utilizing 09-TPS-01 and 
-02 trajectory conditions were analyzed, as these particular trajectories were designed for stressing 
TPS environments [22]. Cases with simulated gap filler protrusions along reference trajectory 08- 
01A were queried for estimates of pressure augmentation factors in the vicinity of step discontinuities 
in the heatshield surface. Comparisons between 05-22 and 09-TPS trajectories are shown against 
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Figure 2.2: CFD Solutions 


the current nominal trajectory 12-GAL-01, where t=0 is defined as the point of entry interface, 
and curves are plotted over the portion of the trajectory where dynamic pressures exceed 850 Pa. 

The CFD database pressures were nondimensionalized using the freestream pressure and dy- 
namic pressure of each solution. The nondimensionalization transforms the forebody pressures into 
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time (s) 


(a) Mach Number 



time (s) 


(b) Dynamic Pressure 


Figure 2.3: Comparison of Mach Number and Dynamic Pressure Along Different 
Trajectories 


pressure coefficients, defined as C v = (p — p 0 o)/Qoo- Modeling the surface pressures in coefficient 
form allows the CFD solutions to be used across a wide range of trajectories with different density 
and velocity profiles. The data were re-interpolated from the 7 block grid to a single-zone grid of 
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clock and cone angles. This alteration simplifies the interpolation of surface conditions without 
the complications of search routines required to handle multiple zones. Each CFD matrix point 
consists of a full surface pressure distribution solution, with 37 clock angles in uniform 5 degree in- 
crements and 61 cone angles with non-uniform increments. These matrix points can be interpolated 
as needed to provide estimates of the pressure distribution at any point on the aeroshell. 

The nominal MEADS pressure model is a tabulation of the nondimensional pressure coefficient 
Cp at each MEADS port mapped over (Af, ctT, <t>) space, where all data are taken from nominal 
Mars entry simulations along the 05-22 trajectory. The dimensionless form is used so that changes 
in absolute surface pressure are scaled according to the freestream dynamic pressure for any desired 
trajectory. The tables were constructed from LAURA results by extracting surface pressure coeffi- 
cients from each (Af,ar) solution at the MEADS port locations, then rotating the solution about 
the body axis by a roll angle increment and repeating as necessary to fill the desired (a,/3) space. As 
MSL’s absolute angle of sideslip is predicted not to exceed 1 degree during the portion of nominal 
flight where freestream dynamic pressure exceeds 850 Pa, the tables have been constructed for each 
(Af, olt) over the range —9 < <f> < 9 in 3 degree increments, providing for more than adequate 
coverage in (3 space. For nonlifting flight after balance mass jettison, the tables are constructed 
over the entire roll angle range from —180 < </> < 180, using the same 3 degree increments in order 
to cover the entire possible range of vehicle attitude. In addition to the surface pressure coefficient, 
the magnitude of its gradient was also calculated using an internal Tecplot function and tabulated 
for use in the pressure uncertainty model, where it is part of a port location uncertainty expression. 


2.3 Experimentation and Simulation 

Note that the pressure model utilized for atmospheric state estimation is based entirely on these 
tabulated CFD solutions. The CFD modeling approach was validated in air in wind tunnel ex- 
periments at Mach numbers of 2.5, 3.5, 4.5, and 10. Residuals between computed and measured 
pressures were within the measurement accuracy obtained in the wind tunnel experiments. 

Wind-tunnel experiments were conducted in the NASA Langley Research Center (LaRC) Uni- 
tary Plan Wind Tunnel (UPWT) with a 6-inch diameter MSL pressure model at Mach numbers 
from 2.5 to 4.5. Two tunnel entries were undertaken, one in 2010 and another in 2011. Test 1835 
was the first UPWT entry and utilized a single model/sting configuration, whereas the second 
entry, test 1855, utilized three model/sting configurations in an attempt to keep the model in the 
same region of the test section regardless of angle of attack. Test 1855 data showed significant 
variability in pressure measurements due to changes in the model location within the test section, 
as shown in Fig. 2.4. Flow surveys show that the UPWT has nonuniform Mach number and flow 
angularity throughout the test section, and that these changes can occur over length scales of the 
same order as the wind tunnel model diameter [23] . A study of pressure error correction techniques 
based on flow survey data were outside the scope of this work. 

UPWT experiments were simulated with the LAURA code, using a calorically perfect gas model 
with the properties of air. The CAD model from which the test article was made was used as the 
surface definition for the computational volume grid, which was composed of 18 blocks and a total 
of 0.77 million points. Cases were run over the angle of attack range 0 < a < 32 in 4 deg increments, 
at nominal Mach 2.5, 3.5, and 4.5 UPWT conditions, all using a constant temperature (300 K) 
viscous wall boundary condition. Grid points were clustered at the bow shock in order to improve 
its placement and reduce its apparent thickness, or “smear,” as shown in Fig. 2.5. 

Experimentation was carried out in the LaRC 31-inch Mach 10 tunnel using the same 6-inch 
diameter model as was used in the UPWT experiments. Experimental results indicated that mea- 
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Figure 2.4: Variability in MEADS Port 1 Measurements due to Changes in Model 
Position 



(a) Unclustered 


(b) Clustered 


Figure 2.5: Grid Lines and Mach Contours at Bow Shock for Unclustered and Clus- 
tered Grids 


sured forebody pressures were extremely sensitive to the model’s location within the test section. 
These data are not used in the pressure or uncertainty models, but are useful to illustrate some of 
the difficulties and challenges associated with pressure model testing in hypersonic ground facili- 
ties, particularly when such experiments are carried out with the intent of validating computations. 
Simulations of the Mach 10 experiments were conducted using the LAURA code, using the same 
grids, models, and assumptions as the UPWT simulations. 
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Chapter 3 


MEADS Transducer Calibration 


A calibration of the MEADS fight system was conducted in order to compute input pressure as a 
function of voltage and sensor /electronics temperature. The MEADS calibration program carries 
two independent calibration efforts. The baseline approach utilizes the industry standard approach 
that has a long heritage and is referred to as the traditional method by the MEADS project. 
This method was used with great success by the SEADS program [5,6] and other FADS programs 
such as the Aeroassist Flight Experiment (AFE) [24]. The second approach is a proof-of-concept 
experiment using a Design of Experiments (DOE) method. This calibration method has never 
before been used to calibrate space fight transducers prior to the MEADS program. The DOE 
approach is documented in [25,26] and is not discussed further in this paper. 

This chapter describes the calibration approach and the analysis of the final calibration data 
using the traditional calibration/analysis approach. The final calibration of the transducers was 
preceded by preliminary thermal testing and acceptance testing. This initial testing led to prelim- 
inary selection and rankings for both Flight and Flight Spare transducer sets. The result of this 
ranking and port assignments are given in Table 3.1. 

3.1 Calibration Methodology: General Overview 

The traditional calibration approach as applied to a FADS type pressure measurement system is 
composed of the collection of pressure transducer performance data and a data analysis process 
originally developed in support of the SEADS experiment. This report briefly describes the calibra- 
tion data collection process and a detailed description of the analysis methodology and assumptions 


Table 3.1: Port Assignments 


Port 

Flight 

Serial Number 

Flight Spare 
Serial Number 

1 

022 

026 

2 

025 

020 

3 

029 

031 

4 

021 

030 

5 

032 

027 

6 

023 

034 

7 

024 

036 
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using final transducer calibration data obtained by the traditional method. 

The traditional method acquires a set of input pressure data vs. transducer output Voltage at 
various temperatures to define transducer sensitivity (scale factor), non-linearity, and bias (zero 
offset) as a function of temperature. The sensitivity and non-linearity data are collected by process 
that uses an increasing/decreasing application of pressure over the operating range of the trans- 
ducer. (In the case of MEADS from zero to 5.5psi.) The number of pressure set points can vary 
depending on project requirements. Typically a total of 11-13 pressure points have proven to be 
sufficient in previous applications. A series of pressure versus Voltage output data sets (curves) 
are obtained at applicable operating temperatures to define the thermal sensitivity of the pressure 
sensitivity and zero offset. The number of temperature set points can also vary. Typically five 
temperature set points have proven to be adequate. 

Other performance characteristics that are determined from the calibration data set are within 
test repeatability and pressure/temperature hysteresis. Repeatability is defined here as the amount 
of change of a measured reading at the same pressure and conditions over a series of pressure 
cycles from zero to full-scale pressure and back to zero again. To ensure no pressure hysteresis is 
introduced into repeatability measurements readings are always taken during an increase in pressure 
or a decrease in pressure but never a mixture of the two in the implementation of the Traditional 
method. 

Typically, due to concerns related to the effect of long space missions where severe environments 
and test requirements can impact transducer performance, multiple full thermal calibrations are 
performed at different stages in the transducer acceptance and test process. These tests then allow 
the quantification of any performance variations that might result from time, environment, testing, 
and/or operation. 


3.2 Implementation of Calibration Approach for MEADS 


The MEADS mission and strain gauge pressure transducers configuration resulted in mission spe- 
cific calibration requirements and a modified Traditional calibration approach. Due to MSL space- 
craft constraints the transducer pressure head containing the strain gauge and the electronics were 
separated and located in regions on the spacecraft with different thermal environments. Addition- 
ally, MEADS required that the transducer performance be maximized at the low-pressure end of 
measurement range. As a result each transducer component - pressure head/strain gauge and elec- 
tronics was calibrated as an independent variable. A thirteen point increasing/decreasing pressure 
sequence with a concentration at the lower end of the pressure scale calibration was performed. The 
final calibration matrix is shown in Table 3.2. Setpoints T2 - T6 are specific to the electronics while 
T7 - T13 are specific to the pressure head strain gauge. Setpoints Til - T 13 were incorporated to 
provide test specific repeatability data due to previous indications relative to transducer thermal 
repeatability and hysteresis. All temperature setpoints make use of the same pressure sequence, 
given in Table 3.3. 

Due to schedule constraints only one set of Traditional calibration data was collected for MEADS 
thus the ability to assess and quantify repeatability as a function of time, environmental effects 
and operation could not be accomplished. 
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Table 3.2: Calibration Temperature Setpoints 


Setpoint 

Number 

Transducer 
Temperature (deg C) 

SSE 

Temperature (deg C) 

T1 

20 

20 

T2 

-80 

60 

T3 

-80 

30 

T4 

-80 

0 

T5 

-80 

-20 

T6 

-80 

-30 

T7 

-50 

0 

T8 

-80 

0 

T9 

-100 

0 

T10 

-125 

0 

Til 

-100 

0 

T12 

-80 

0 

T13 

-50 

0 

T14 

20 

20 


Table 3.3: Pressure Sequence 


Sequence Number 

Pressure (psi) 

1 

0.00 

2 

0.25 

3 

0.50 

4 

1.00 

5 

2.00 

6 

3.30 

7 

4.50 

8 

5.50 

9 

3.90 

10 

2.50 

11 

1.50 

12 

0.75 

13 

0.00 
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3.3 Calibration Analysis: Methodology and Assumptions 


The traditional calibration method uses a least-squares fit for each transducer to a model of the 
form [6] 

p = a 0 (Ti,T 2 ) + ai(Ti,T 2 )l/ + a 2 (T 1 ,T 2 )y 2 (3.1) 


where T\ is the pressure transducer temperature, T 2 is the SSE temperature, V is the Voltage 
output, and p is the pressure. The coefficients ao, aq, and a 2 are termed the bias, sensitivity, and 
non-linearity, respectively. These coefficients are determined from a least-squares fit of pressure 
vs. Voltage at specific temperature setpoints using a simple increasing/decreasing applied pressure 
to create the Voltage response. By varying the temperature setpoints, a database of calibration 
coefficients can be constructed, and an empirical model developed that uses spline interpolation to 
look up the values of these coefficients and functions of the two temperatures. 

The specific calculations performed at each temperature setpoint is provided in the following. 
First, the matrix F is defined as 

r i ^ v? i 


F = 


1 V n 



(3.2) 


where n is the total number of pressure points at each temperature setpoint. Similarly, the array 
of pressure measurements is written as 


P = 


Pi ' Pn 


1 T 


(3.3) 


The coefficient array is written as $, = 
coefficient array is given by 

s = 


1 T 


a o a i a 2 
-l 




F t p 


The least-squares estimate of the 


(3.4) 


An important byproduct of the least-squares solution is the covariance matrix of the coefficient 
estimate, given by 


Co v U = a 2 F T F 


-i 


(3.5) 


where a c is the Mensor pressure controller accuracy. Along with the values of the calibration 
coefficient, the uncertainties in the coefficients are also stored as a function of the two temperatures. 
These uncertainties can be used to build a calibration error model to assess system performance. 
The total system uncertainty model from the calibration also needs to consider the temperature 
uncertainties and hysteresis. To this end, define 


z = 


T\ T 2 V a o oi a 2 


1 T 


(3.6) 


It follows 
ance is given 


from a linear covariance transformation that the total pressure prediction error vari- 

by 






T 

+ d lo 


+ '<V 2 + al 2 V 4 


(3.7) 


where 


n z 


< 7 ^ 0 0 0 

0 <jy 2 0 0 

0 0 <4 0 

0 0 0 Cov (x) 


(3.8) 
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The quantities a^ 0 , , and a^ 0 account for non-repeatability and hysteresis in the bias, sensi- 

tivity, and nonlinearity coefficients in the calibration equation, respectively. These contributions to 
the uncertainty model are computed from repeated points in the calibration temperature setpoints. 
In summary, the uncertainty predictions include non-repeatability/hysteresis, calibration error, 
temperature error, and Voltage error/noise. Note that in this model the temperature uncertainty 
due to thermocouple bias and noise are lumped together into the total temperature uncertainty. 

The in-flight zero obtained prior to Mars entry has the effect of significantly reducing the 
uncertainties in the pressure data. Essentially, the in-flight zero nearly removes the uncertainties 
due to bias, but the slope and nonlinearity uncertainties remain, and a small residual bias error 
remains due to the slope of the bias curve as a function of temperature, coupled with temperature 
uncertainty and noise. An example of the effect of the in-flight zero is shown in Figure 3.1 based on 
data from Transducer 022 at temperature setpoint T4 (Transducer temperature -82.9 C and SSE 
temperature 9.7 C). Figure 3.1(a) shows the pressure vs. Voltage curve from the calibration data 
for the (true) Mensor controller pressure and the estimates pressure with and without the in-flight 
zero applied. On this scale, the transducer behaves very linearly with a small correction / bias shift 
due to the in-flight zero. Figure 3.1(b) shows the pressure error with and without the correction. 
On this scale, a hysteresis loop and the bias correction due to the in-flight zero can be observed. 

To develop an uncertainty model, let Tf and T| be the transducer and SSE temperatures, 
respectively, and V* be the measured Voltage at the time of the in-flight zero prior to atmospheric 
entry. Then the pressure error is 

5p = a 0 (7?, T|) + 0l (I?, T 2 *) V* + a 2 (7?, T 2 *) (V*) 2 (3.9) 

It follows that the adjusted pressure is given by 

p = p — dp (3.10) 


By substituting Eq. (3.9) into Eq. (3.10) and linearizing for small temperature errors, the 
pressure measurement uncertainty model following the in-flight zero is given by 



where a ni , ai n , cr n2 , and cq, 2 are the thermocouple noise and bias uncertainties for temperature 
measurements T\ and T 2 , respectively. 


3.4 Analysis of Final Calibration Data 

Data from the final calibration of the MEADS/SSE system were analyzed according to the method- 
ology described in the preceding section. The traditional data points collected during the calibration 
are shown in Table 3.2. The desired vs. actual temperature setpoints acquired in the calibration 
are shown in Fig. 3.2 for the Flight and Flight Spare calibrations. Additional data were collected 
for the DOE calibration, which included several “pressure profiles” that attempt to simulate the 
pressures during Mars entry. No data were excluded from these datasets, although the ambient test 
points were analyzed but are not incorporated into the flight data reduction database since these 
points are not within the flight temperature space. No other data was utilized for the purposes of 
this calibration report, including ambient performance tests and acceptance testing data. 
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-2 0 2 4 6 8 10 12 

Voltage [mV] 

(b) Pressure Error 

Figure 3.1: Effect of In-Flight Zero 


Recall that replications are made involving temperature setpoints (T1,T14) at ambient, (T4, 
T8, T12) at -80/0, (T7, T13) at -50/0 and (T9, Til) at -100/0. These replications allow the 
repeatability and hysteresis of the transducer to be evaluated. Both the Flight and Flight Spare 
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transducers make use of the same temperature set points in order to enable direct comparisons 
between both systems. 

The non-repeatability/hysteresis uncertainties were calculated from the repeated temperature 
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Table 3.4: Flight Transducer Non-Repeatability/Hysteresis Uncertainties 


Transducer 

Serial 

Number 

do Bias 
Uncertainty 
[%FS] “ 

(j i Sensitivity 
Uncertainty 
[% FS/rnV] 

a 2 Nonlinearity 
Uncertainty 
[% FS/rnV 2 ] 

Within Temp. 
Repeatability 
[% FS] 

Between Temp. 
Repeatability 
[% FS] 

020 

0.348 

7.98e-3 

4.12e-4 

0.051 

0.347 

021 

0.838 

1.06e-2 

6.126-4 

0.060 

0.735 

022 

1.233 

7.38e-3 

5.816-4 

0.060 

1.078 

023 

0.658 

5.04e-3 

4.916-4 

0.067 

0.573 

024 

0.470 

3.62e-3 

3.03e-4 

0.042 

0.414 

025 

0.769 

8.87e-3 

4.03e-4 

0.065 

0.684 

026 

0.407 

7.50e-3 

4.16e-4 

0.048 

0.407 

027 

0.225 

9.31e-3 

6.11e-4 

0.042 

0.225 

029 

0.914 

5.79e-3 

3.50e-4 

0.059 

0.801 

030 

0.633 

8.23e-3 

4.83e-4 

0.071 

0.613 

031 

0.428 

8.48e-3 

6.11e-4 

0.064 

0.416 

032 

0.619 

9.96e-3 

6.75e-4 

0.054 

0.545 

034 

0.937 

3.35e-2 

3. Ole-3 

0.046 

0.960 

036 

0.945 

1.25e-2 

9.25e-4 

0.040 

0.956 


setpoints listed above. The results for the transducers are shown in Table 3.4. This table shows 
the 1-ex uncertainties in the calibration equation coefficients used in the uncertainty modeling, 
with pressure units in percent full scale. It is possible to generate rankings based on any of these 
uncertainties, but it is not recommended, since the total system uncertainty is a preferable metric 
as it includes all uncertainty sources together. However, several interesting observations can be 
made about this data. Namely, transducer 024 has the lowest uncertainty in all three coefficients 
among the Flight set. Many transducers from the Flight Spares set (020, 026, 027, and 031) have 
slightly better bias uncertainty than 024, but 024 exhibits the best sensitivity and nonlinearity 
uncertainty among all transducers. Transducer 034 is an outlier in the sense that its sensitivity 
and nonlinearity uncertainties are an order of magnitude larger than all other transducers. Also 
shown in these tables are the Within and Between Temperature Repeatability uncertainties. These 
uncertainties are utilized by the DOE method and are shown here simply for comparison purposes. 
Details on how these uncertainties are calculated can be found in [27]. Interestingly, the between 
temperature repeatability trends correlate with the bias uncertainties. 

The quadratic calibration equation coefficients were calculated from the least-squares fit to the 
pressure/ Voltage data across all temperature setpoints in the region of interest to flight data re- 
duction (i.e. T2-T13). Next, the repeated points were averaged together and extrapolations were 
conducted, treating the transducer temperature and SSE temperature as independent variables, to 
create data tables suitable for 2-Dimensional look-up routines to determine the pressure as a func- 
tion of Voltage, transducer temperature, and SSE temperature. The look-up tables are tabulated 
in Appendix A. For the results described in this document, the Matlab-based griddata routine was 
utilized to perform the look-ups, using the cubic interpolation method. The linearizations required 
for calculating the uncertainty estimates are computed using a centered divided difference method 
with 0.5 deg C temperature perturbations. 

The resulting calibration model was utilized to predict the DOE and Pressure Profile data (with 
and without the simulated in-flight zero) to asses performance of the model itself as well as the 
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Table 3.5: Flight Transducer Non-Repeatability/Hysteresis Uncertainties 


Transducer 

Serial 

Number 

DOE 
Actual 
RMS Error 
[% FS] 

DOE 
Predicted 
RMS Error 
[% FS] 

Profile 
Actual 
RMS Error 
[% FS] 

Profile 
Predicted 
RMS Error 
[% FS] 

Profile 
Actual 
RMS Error 
[% FS] 

Profile 
Predicted 
RMS Error 
[% FS] 

020 

0.385 

0.408 

0.217 

0.404 

0.214 

0.080 

021 

0.928 

0.805 

0.841 

0.801 

0.229 

0.074 

022 

1.326 

1.200 

1.316 

1.197 

0.245 

0.067 

023 

0.660 

0.655 

0.600 

0.654 

0.138 

0.085 

024 

0.460 

0.459 

0.392 

0.458 

0.110 

0.045 

025 

0.816 

0.756 

0.750 

0.754 

0.175 

0.045 

026 

0.705 

0.497 

0.174 

0.497 

0.182 

0.066 

027 

0.251 

0.266 

0.252 

0.263 

0.194 

0.087 

029 

0.962 

0.891 

1.008 

0.890 

0.173 

0.079 

030 

0.630 

0.708 

0.414 

0.707 

0.302 

0.042 

031 

0.356 

0.470 

0.279 

0.465 

0.196 

0.076 

032 

0.680 

0.610 

0.677 

0.606 

0.171 

0.049 

034 

1.277 

1.055 

0.704 

1.038 

0.667 

0.102 

036 

0.884 

1.039 

0.624 

1.037 

0.467 

0.057 


uncertainty predictions based on the linear covariance transformations. The results are tabulated 
in Table 3.5. These results show that the predicted uncertainties agree with the actual uncertainties 
reasonably well. The in-flight zero uncertainty predictions do not match as well as the other cases, 
but this difference could be due to how the profiles were executed during the calibration. Further 
analysis is required to show this is the case. 

Based on these results it is possible to generate rankings of transducer performances based on 
the in-flight zero results, which are the most representative of actual system performance. Sorted 
in ascending order, the recommended flight system ranking is 024, 023, 032, 029, 025, 021, and 
lastly 022. Since the best two transducers are currently assigned to ports 6 and 7, recommendations 
were made to re-order the port assignments in order to improve the dynamic pressure and angle of 
attack performance. Due to schedule constraints it was not possible to do so. 


3.5 Uncertainty Analysis 

An analysis of uncertainties from the final calibration data was conducted to examine transducer 
performance relative to system requirements of 0.5% accuracy of pressure reading. The analysis 
described in this section is based on the 09-HOL-01 nominal trajectory. The nominal port pres- 
sures along this trajectory are shown in Figure 3.3. The 850 Pa dynamic pressure boundaries are 
indicated, as well as the 5 psi full scale pressure threshold. 

The calibration uncertainties along the reference trajectory are shown in Figure 3.4. Note that 
the error uncertainties are shown as a percent of reading. These uncertainties include the transducer 
calibration uncertainties alone, with all other error sources ignored and the in-flight zero applied. 
The errors are based on the 09-HOL-01 nominal trajectory. These curves have been generated 
using the expression for the error variance given in Section 3.3 with appropriate terms set to zero. 


22 



40 


35 i 


30 

'(7 

CL 

i 25 

a> 

&. 

3 

w 20 

w 

o 

L- 

Q. 

r is 

o 

a 

10 


0 L 

0 


Port 1 
Port 2 


Port 3 
Port 4 



50 100 150 200 250 

Time Since Entry [s] 


Figure 3.3: 09-HOL-01 Nominal Pressures 


The reduced expression is 

4 = (< + ~<) V 2 + (4 + 4 ) ^ ( 3 - 12 ) 

The error transformed into percent of reading is simply 

%Error = ^ x 100 (3.13) 

The uncertainty reduces to a constant value at zero pressure due to the limit being a ratio of 
two slopes. To a first-order approximation (equating 10 mV to 5 psi), an expression for the error 
as a function of pressure is 

Ap = eip + £2P 2 (3-14) 

where ei = 2Aai and C 2 = 4Aa2- The percent error reduces to 


%Error 


^ ) x 100 = (ei + e 2 ) x 100 = 100ci + 100e 2 p 
P ) 


(3.15) 


This result indicates that the percent error is a linear function of pressure, reaching a constant 
value at zero pressure and maximum value at peak pressure. Note that the percent error results 
are less than 0.5% for all transducers except 034 and 036. Transducer 036 has less than 0.5% error 
over most of the trajectory but exceeds slightly during the peak dynamic pressure regime. All other 
transducers meet the requirement of less than 0.5% error over the entire trajectory. 
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Chapter 4 

Uncertainty Modeling 


At the MEDLI project outset, no top-down error budget was formulated to determine the hardware 
requirements. To meet the available schedule at the time, analysis or engineering judgment was 
used to identify the sources of error that were both (1) expected to contribute significantly to 
errors in the returned data, and (2) controllable through manufacturing. The MEADS Level 4 
requirements were written to minimize uncertainties in the final data return. The hardware was 
then produced, calibrated, and integrated to the Level 4 requirements. 

In parallel with MEADS hardware development, an error allocation was devised for the angle of 
attack and angle of sideslip data products. This error budget hierarchy is shown in Figure 4.1. The 
Level 2 requirement for angle of attack knowledge (PS-372) is +/- 0.5 degrees 3-sigma, shown in the 
top box of the hierarchy. The contributors to the uncertainty come from the transducer assembly, 
the Sensor Support Electronics, the Aeroshell manufacturing, and the system time latency. Each 
box lists the uncertainty requirement on that specific contributor, or the as-built value achieved, 
if applicable, and the error in angle of attack that results from that uncertainty. The following 
sections explain the analysis that produced these results. 

The hierarchy for the angle of sideslip is identical except for the pressure transducer contribution 
to the errors; each of the seven transducers is unique, and sideslip takes advantage of the two ports 
and transducers that are symmetric about the vehicle centerline. 

The impact of individual error sources and the effects of selected failure modes are examined 
in this chapter. These error sources include port location knowledge, sensor noise and quantiza- 
tion, telemetry time stamp errors, pneumatic lag, thermal transpiration, and system calibration 
uncertainties. Errors due to CFD pressure distribution uncertainties are also described in this chap- 
ter. The associated error models are reviewed in the following sections before showing sensitivity 
simulation results and the overall error budget results compared to the project requirements. 
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Angle of Attack 
Required: +/- 0.5 deg 3sigma 
RSS of measured: 0.38 deg 


Transducer Error 
Required: +/- 0.5% of reading 


f— 

SSE Error 
Requirement: None 


Manufacture/Geometry Error 
Alpha error: 0.05 deg 


Time Latency Error 
As-Built: +/- 25us 3a 

Alpha error: 0 345 deg 



Alpha error: 0.14 deg 


Deformation not considered 


Alpha error: 0 deg 


Calibration Error (w/ IFZ) 
Alpha error: 0.33 deg 


Transducer Temperature Error 
Alpha error: 0.10 deg 


Thermal Transpiration Error 
Required: +/- 0.5% of reading 
Alpha error: ~0.0 deg 


Pneumatic Lag Error 
Path length + Path diameter 
Alpha error: < 0.01 deg 


SSE Quantization 
Alpha error: ~0.0 deg 


SSE Noise 
Alpha error: 0.14 deg 


SSE Temp Uncertainty 
Alpha error: 0.02 deg 


Port Location Knowledte 
Design: 0.05in. 3sigma 
Alpha error: 0.05 deg 

Aeroshell Deformation 
Requirement: < 0.5 in. 
Achieved: 0.3 in. 

Alpha error: Not considered 


1 553 Bus Latency Uncrt. 
Design: +/- 21.5ps 3a 






SSE Flight Oscillator Variation 
As built: +/- 8pS 3a 




r 


RCE+DPAM Flight 
Oscillator Variation 
As built: +/- 8us 3o w 


_y 


r- 

f ! + [ ] + Used in uncertainty analysis 

i— Not included in analysis or not 

L presented (contribution negligible) 

Not included (out of scope for MEADS SAR) 


Figure 4.1: Angle of Attack Error Budget Hierarchy 


4.1 Hardware System Error Models 

4.1.1 Port Location Uncertainties 

Errors in the knowledge of the port location after installation are cast as random longitudinal 
perturbations at constant clock angle and a random tangential perturbation at constant cone angle. 
The errors are specihed as uncorrelated zero- mean Gaussian random numbers. The measurement 
system is accurate to within 0.04. 


4.1.2 Sensor Noise and Quantization 

Sensor noise is modeled using the equation 

3p = iPmax Pmin ) Vp ( 4 . 1 ) 

where pmax and p m in are the maximum and minimum pressure range of the sensor system, respec- 
tively. The variable r] p is a uniform random variable between ±N/2 n where N is the number of 
counts of noise and n is the system effective bit size. 
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Sensor quantization can be modeled using the equation 


( Pmax Pmin , T , 

p = lilt 

\ 2 n 


P 

L \Pmax Pmin 


Finally, sensor saturation can be modeled as p = max [min ( p,p m in ) ,Pmax}- 


(4.2) 


4.1.3 Time Tag Error 

Relative errors in the telemetry time stamp between measurements at different port locations 
introduces a distortion effect that corrupts the pressures used as input to the basic data reduction 
algorithm. Timing errors include both random and systematic components. The systematic timing 
errors can the adjusted in pre-processing of the flight data, leaving only the random error to be 
considered in error budgeting. These errors are simulated by first casting a uniform random number 
for each port in the range of ±25 ps, and then shifting the truth state used to calculate the measured 
surface pressures. 


4.1.4 Pressure Leak Rate Error 

The MEADS pressure transducers cannot be tested with applied pressure once they are mounted on 
the MSL aeroshell, since contact with the PICA thermal protection system material is prohibited. 
Therefore, the project developed and tested an installation process that minimizes leaks along the 
pressure path. Pressure leak analysis is described in 28,29. 

Pressure system leaks can be modeled as a multiplicative loss on the measured pressure at the 
transducer relative to the port pressure. Since pressure leaks are a one-sided distribution, the errors 
can be reformulated as a two-sided distribution with half the total error by adjusting the pressure 
up, assuming the mean value of the leak rate. The adjusted pressure measurement is given by 

P=(! + ^P)p (4.3) 


where L max is the maximum expected leak rate, assuming that the leak rate is uniformly distributed. 
Following this adjustment, the pressure leak errors are uniformly distributed in the interval 


£ eU 


Lmax 

2 


Lmax 

2 


(4.4) 


It then follows that the simulated pressure errors due to leaks are given by 


p = (1 + £)p 


(4.5) 


4.1.5 Vibration Error 

Tests were performed to assess pressure error due to vibration during EDL. The tests were performed 
with 4 MEADS transducers (S/N 002, 028, 033, and 035) at 1.33 grms vibration conditions during 
a simulated pressure profile. Three vibration tests were performed for each axis of the transducer, 
where the y-axis test is directed along the pressure path. The measured transducer outputs from 
these tests are shown in Figures 4. 2-4. 4. 

The noise due to vibration was then estimated from this test data. First, an “in-flight” zero was 
applied to the data to remove bias, and then residuals were computed between the true pressure 
and the measured pressure. These residuals were filtered using 1 Hz optimal Fourier filter. The 
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Figure 4.2: X-Axis Vibration Test Data 


Table 4.1: Vibration Errors, la % Full Scale 


Transducer 

Axis 

Stationary Noise 

Unfiltered Error 

Filtered Error 

PI 

X 

0.0135 

2.5576 

1.245 

P2 

X 

0.0097 

2.7153 

1.2752 

P3 

X 

0.0124 

3.8236 

1.7622 

P4 

X 

0.0104 

3.3704 

1.6166 

PI 

Y 

0.0102 

3.8372 

1.8969 

P2 

Y 

0.0072 

10.029 

4.4757 

P3 

Y 

0.0089 

3.3062 

1.559 

P4 

Y 

0.0076 

5.2139 

2.5568 

PI 

Z 

0.0167 

0.4979 

0.3481 

P2 

Z 

0.0346 

0.6368 

0.3618 

P3 

Z 

0.0143 

0.4401 

0.3044 

P4 

Z 

0.0174 

0.4994 

0.3291 


errors due to vibration were then computed from these filtered residuals. Error statistics from these 
residuals are shown in Table 4.1. 

From these results, an uncertainty model of the form 

(y noise = ( 4 - 6 ) 

H.max 

is introduced, where cr v ib e is the total measured error, <tq is the stationary noise uncertainty, and 
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Time [s] 


Figure 4.3: Y-Axis Vibration Test Data 


@ noise is the noise uncertainty due solely to vibration. The vibrational noise term is then scaled by 
the ratio of instantaneous dynamic pressure to the maximum dynamic pressure. 

Next, a study of this noise model for different filtering frequencies was conducted. First, the 
errors introduced by filtering perfect data was studied. A 1 Hz filter and a 0.1 Hz filter were 
applied to the data. In this case, some true high frequency vehicle dynamics are removed by low 
pass filtering. This behavior is illustrated for the Port 5 pressure in Figure 4.5. Figure 4.5 shows 
the entire pressure profile, Figure 4.6 shows a detail view of the pressure at maximum dynamic 
pressure. Here it can be seen that the 0.1 Hz filter effectively removed all rigid body dynamics of 
the vehicle. Figure 4.7 shows the errors in pressure that arise due to filtering the truth data. 

Figures 4.8 and 4.9 show a similar set of results for the case with 0.5% (of full scale) measure- 
ment noise at peak dynamic pressure. Figure4.10 shows the errors between unfiltered and filtered 
pressures at 1 and 0.1 Hz. Here, the 0.1 Hz filter has large errors because it has filtered the rigid 
body dynamics. The 1 Hz filter produces lower errors because it filters some of the noise while 
retaining the rigid body motion of the capsule oscillations. 

A more stressful case is shown in Figure 4.11-4.13, in which the vibrational noise is set to 5% of 
the full scale pressure at peak dynamic pressure. In this case, the noise is large enough in amplitude 
that it dominates over the capsule’s rigid body motions. As a result, the 0.1 Hz filter produces 
a better estimate of the true pressure. The results of these three experiments are summarized in 
Table 4.2, which show the 1 -a pressure errors over the entire trajectory for each case. These results 
indicate that the choice of filter cutoff frequency should be guided by the expected vibration noise 
environment. This is further illustrated in Figures 4.14-??, which shows power spectral densities 
for each of the three test cases, and for unfiltered, 1 Hz filter, and 0.1 Hz filter data. From these 
plots, it can clearly be seen that the vibration noise dominates over the rigid body modes as the 
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Figure 4.4: Z-Axis Vibration Test Data 


noise level increases. Given these results and the expected noise environment based on vibration 
tests, the 1 Hz filter is recommended for use on MEADS data reduction. 

Table 4.2: Pressure Errors over Trajectory, la % Full Scale 



No Noise 

0.5% FS Noise 

5.0% FS Noise 

Unfiltered 

0.0 

0.1596 

1.5817 

1 Hz filter 

0.0471 

0.0932 

0.7289 

0.1 Hz filter 

0.1172 

0.1193 

0.2026 
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Port 5 Pressure [psi] Port 5 Pressure [psi] 



Figure 4.5: Port 5 Pressure: No Noise 



Time [s] 


Figure 4.6: Port 5 Pressure: No Noise (Detail) 
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Port 5 Pressure [psi] Port 5 Pressure [psi] 



Time [s] 


Figure 4.7: Pressure Error Due to Filtering 



Figure 4.8: Port 5 Pressure: 0.5% Full Scale Noise at Maximum Dynamic Pressure 
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Figure 4.9: Port 5 Pressure: 0.5% Full Scale Noise at Maximum Dynamic Pressure 
(Detail) 
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Figure 4.10: Pressure Error Due to Filtering 
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Figure 4.11: Port 5 Pressure: 5% Full Scale Noise at Maximum Dynamic Pressure 



Figure 4.12: Port 5 Pressure: 5% Full Scale Noise at Maximum Dynamic Pressure 
(Detail) 
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Figure 4.13: Pressure Error Due to Filtering 
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Figure 4.14: Port 5 Power Spectral Density: No Noise 
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Port 5 Spectral Density [dB] Port 5 Spectral Density [dB] 
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Figure 4.15: Port 5 Power Spectral Density: 0.5% Noise 
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Figure 4.16: Port 5 Power Spectral Density: 5% Noise 
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4.1.6 Pneumatic Lag 

The MEADS system response model, or lag model, captures the dynamics of the pressure mea- 
surement from the surface port location to the pressure transducer (the sensor). In general, these 
characteristics are a function of the line properties (tube diameter and length), fluid properties 
(density, dynamic viscosity, etc.), as well as flight conditions (pressure and temperature at the 
port and temperature at the transducer). The system response models are used to accomplish 
two objectives: firstly, to accurately model the response so that pre-flight simulation can be used 
to assess hardware, algorithms, and telemetry rates, and secondly, to characterize the system re- 
sponse so that measured pressures can be corrected to actual surface pressures for use in generating 
atmospheric data estimates. 

High-fidelity response models based on the Navier-Stokes equations are available in the liter- 
ature, in particular the infinite-order frequency-domain Bergh-Tij deman response model [30] and 
variations on this model for branched systems developed by Whitmore [31]. In these models, the 
pneumatic system is broken up into a series of segments, each with constant line and fluid prop- 
erties. The end to end frequency response from port to transducer becomes the product of the 
individual response functions for each segment. Entrapped volume in any of the segments can also 
be captured. In this manner changes in tubing diameter as well as temperature gradients can be 
explicitly modeled. The response is evaluated across a range of input frequencies and the result is 
a complex frequency response that can be represented as magnitude and phase. 

The frequency domain results for each flight condition are not readily suited either for time 
domain simulation or for direct analysis of time histories of measured pressure. For these purposes, 
a low-order response model must be utilized. For simple 2nd order systems, this response can be 
characterized by a natural frequency and damping, with transfer functions of the form 


Pq S 2 + 2 £ui n + U! 2 

where Pq is the applied pressure, Pl is the lagged pressure, and ( and u n are the damping ratio 
and natural frequency, respectively, and s is the Laplace variable. 

A variety of 2nd order models based on simplified fluid mechanics relations are available in the 
literature with varying assumptions and classifications [32] . A new 2nd order response model based 
on first principles is presented in [32], This model is derived directly from the Bergh-Tij deman 
equations using order reduction techniques and was found to be more accurate than other 2 nd order 
models. The natural frequency and damping coefficients of this model are 


(4.7) 


u 


2 

n 


c 


A c a 2 /Lv 

(■ LA c /2v ) [l + ( 1 / 12 ) ( 7 LR d ) 2 / (a s po) 2 


+ 1 + (1/6) {'yLR d ) 2 / (a s po) 2 


u n 

2 


7 Pd 
Po 


vL 1 
Acol + 2 


(4.8) 

(4.9) 


where L is the tube length, A c is the cross-sectional area, v is the transducer entrapped volume, 
po is the fluid density in the tube, c is the speed of sound, and 7 is the ratio of specific heats. The 
parameter R ^ is known as the acoustic impedance and can be calculated as a function of Reynolds 
number to capture laminar, transitional, or turbulent flow conditions, as described in [32]. 

Another approach to time domain analysis is look for a low-order system that matches the model 
response. This assumed form of this equivalent system is chosen to capture the modeled dynamics 
and to facilitate efficient simulation. Numerical optimization techniques have been utilized to fit 
transfer function models to the high-fidelity frequency response model, including up to 4th order 
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systems with pure time delay. These results were not found to significantly improve the 2nd order 
model and thus were not used to model lags for MEADS error analysis. 

Due to computational complexities of simulating the tube and transducer response, these models 
were employed in a stand-alone manner in order to assess impacts to system performance. These 
studies include both the high-fidelity frequency domain models suitable for analyzing specific flight 
conditions as well as 2nd-order approximate models suitable for analyzing time histories. 

The first study is an assessment of the port diameter of 0.1 in. Three flight conditions were 
explored using the high-fidelity frequency response model, including a high Mach/low pressure 
case, the maximum pressure case, and a low Mach/low pressure case. The frequency response for 
these three conditions are shown in Fig. 4.17 and 4.18. The frequency response of the 2nd order 
model is also shown for comparison. Table 4.3 shows several important metrics of the frequency 
response for these three flight conditions. Here, the time constant and the time delay (computed 
from phase loss at 1 Hz frequency) are shown corresponding to the high-fidelity (HF ) model as well 
as the 2nd order model. While differences are apparent between the HF and 2nd order models, the 
trends and orders of magnitude are similar. These results agree to within a factor of 2 or better with 
an independent implementation using a reduced-order Bergh-Tij deman model described in [33]. 
These results indicate that the response time is acceptable even for the worst-case lag time at the 
low Mach, low pressure condition. 



A time history of the pressure difference computed using the 2nd order model is shown in 
Fig. 4.19. If the input signal is of sufficiently low frequency content, then these lag effects are 
completely deterministic in nature and can effectively be corrected for by inverting the transfer 
function to solve for the pressure signal at the port. Assuming no model uncertainties, then the 
lag effects can be corrected without error. However, errors can arise due to the presence of tube 
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Frequency [Hz] 


Figure 4.18: Tube Response: Phase 


Table 4.3: Response Model Conditions and Results 


Condition 

Moo 

goo, Pa 

P, Pa 

t, sec (HF) 

r, sec (2nd) 

At, sec (HF) 

At, sec (2nd) 

1 (High Mach) 

29.0 

850.0 

1630.0 

0.34-10 -3 

0.72-10 -3 

0.70-10 -3 

0.75-10 -3 

2 (Max ^oo) 

16.8 

14015.3 

26788.5 

0.45T0 -3 

0.33T0 -3 

0.55-10 -3 

0.14-10 -3 

3 (Low Mach) 

3.1 

850.0 

1595.7 

1.31-10 -3 

0.78T0 -3 

2.01-10 -3 

1.09-10 -3 


geometry uncertainties. Therefore it is of interest to understand the sensitivities to diameter and 
length uncertainties. 


Table 4.4: Response Model Port Diameter Results 


Port Diameter, in 

r, sec (HF) 

r, sec (2nd) 

At, sec (HF) 

At, sec (2nd) 

0.07 

6.13-10 -3 

0.86-10 -3 

6.61 ■ ur 3 

1.21-10~ 3 

0.04 

5.69-10 -2 

0.96-10 -3 

5.49T0 -2 

1.36-10” 3 

0.02 

2.62 

1.04-10 -3 

0.22 

1.47-10” 3 


Table 4.4 shows the results of a study varying the port diameter for the worst-case flight 
condition (low pressure, low Mach) discussed above. Here, the port diameter is varied to as low 
as 0.02 in. These results indicate that the 0.07 in case may be acceptable but the 0.04 in case 
and below quickly become unacceptable. Figures 4.20 and 4.21 show the frequency response of the 
high-fidelity and 2nd-order models for these port diameters. Discrepancies between the HF and 2nd 
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Figure 4.19: Nominal Pressure Lab 
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order models tend to increase with decreasing port diameter. Arc jet tests conducted for various 
port diameters indicate that the 0.10 in case is acceptable, therefore the diameter has been kept as 
originally designed. However, these results serve to indicate the sensitivity to the port diameter, 
therefore the diameter uncertainties should be kept as small as possible. The port 5 (worst-case 
port) pressure lag error due to diameter uncertainty is shown in Fig. 4.22 for a range of diameter 
uncertainty up to 0.03 in. In all cases the lag error is small (note that the errors are shown in % of 
reading) . 


Table 4.5: Response Model Tube Length Results 


Tube Length, in 

r, sec (HF) 

r, sec (2nd) 

At, sec (HF) 

At, sec (2nd) 

+0.1 

1.22- 1 0 -3 

0.80T0 -3 

2.03T0 -3 

LIMO -3 

+0.25 

1.26-10" 3 

0.82-10 -3 

2.08T0 -3 

1.14-10 -3 

+0.5 

1.31T0 -3 

0.84-10" 3 

2.15T0 -3 

1.18T0 -3 


Next, a study was conducted to assess the sensitivity to tube length. Three cases were in- 
vestigated using at the low pressure, low Mach flight condition. Table 4.5 summarizes the effects 
perturbing the total tube length by increments of 0.1 in, 0.25 in, and 0.5 in. In each case the 
time constant and time delay metrics do not change considerably. Figures 4.23 and 4.24 show the 
frequency response results for these perturbations. These results indicate that the pressure system 
is not highly sensitive to the tube length. Figure 4.25 shows the time history of the lag error for 
port 5 induced by length uncertainties of up to 0.5 in. As in the previous case, these errors are 
acceptably small. 
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Figure 4.22: Port 5 Lag Sensitivity 


Based on the system requirements of 1.41 mm length uncertainty and 0.254 mm diameter 
uncertainty, a worst-case lag error over the entire trajectory was determined from these sensitivity 
results to be 6.94e-04% of reading at each port, in a 3 a sense. 
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Figure 4.23: Length Sensitivity: Magnitude 
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Figure 4.25: Port 5 Lag Sensitivity 


4.1.7 Thermal Transpiration 


Under rarefied flow conditions, transpiration induced effects within a tubing system can cause the 
pressure measured by the transducers to be significantly different from the surface pressure. The 
effects can be divided into those caused by temperature gradients between the transducer and the 
surface, and those caused by heat transfer and/or shear stress. 

A study of the available literature yielded three models to account for transpiration effects; two 
models for simulating the difference between measured and orifice pressure due to thermomolecular 
flow effects and one model for simulating the difference between the surface and orifice pressure due 
to shear stress and/or heat transfer. The first transpiration model is provided in [34], which is an 
application of thermomolecular flow theory developed by Knudsen in [35,36]. Knudsen’s general 
formula for the thermomolecular pressure effect is 


dp 

dT 


1 P 

(1 + b'r/X) 2 2T’ 


where b' 


b 0 


1 + b\r / A 
1 + b 2 r/X 


(4.10) 


where r is a physical length scale, A is the mean free path, P is the pressure, T is the temperature, 
and bo, b \ , and b 2 are empirical constants. The empirical constants have been obtained by Knudsen, 
using least-squares fit to data acquired with hydrogen in glass tubes [35]. The results are bo = 2.46, 
b\ = 3.15, and b 2 = 24.6. 

Assuming that the temperature in the tubing between the surface and the transducer is known 
as a function of length x, Eq. (4.10) can be written as 


dp ( x ) p (x) 

dT (x)' 

« (p, T, D f + b 0 K (p, T, D) 

dx 2 T (x) 

dx 

_k(p, T, R ) 2 + (b 0 + bi) k (p, T,D) + b 1 b 2 _ 
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Table 4.6: Gas Constants for Carbon Dioxide 


Constant 

Value 

Rg, J/kg-K 

191.8 

Mo , Pa-s 

14.8T0 -6 

T 0 , K 

293.15 

C, K 

240.0 


where the parameter k is the Knudsen number, given by 



irRgT 

2 


(4.12) 


where R g is the specific gas constant. The viscosity p can be calculated as a function of temperature 
using Sutherland’s relation [37] 


M = Mo 


/7o + C\ 

\T + C ) 


{ T \ 


(4.13) 


The constants associated with carbon dioxide are summarized in Table 4.6. 

The second transpiration model is based on an algorithm described by Potter. [38]. This model 
is based on the theoretical formulation of [39] and empirical data obtained in [40]. The empirical 
data were obtained in a sequence of experiments involving a wide range of tube lengths, diameters, 
orifice shapes, and gases, including hydrogen, helium, nitrogen, and argon. The model is 


dp ( x ) 
dx 


p (x) 
2 T (x) 


dT (x) 
dx 


i + 00 + 01 

n (p, T, D) K (p, T, D) + c 2 /k (p, T, D ) J 


-2 


(4.14) 


where Co = 0.275, c\ = 0.625, and c 2 = 24.0. 

Note that both the Knudsen and Potter transpiration models predict that ( dp/dx ) — > 0 as 

n — y 0 . 

The last model investigated for this paper is an augmentation of the Potter model [38] that 
includes non-equilibrium effects due to heat transfer and shear stress. This model is described 
in [41]. The solution relies on empirical fits to data obtained in a wide range of experiments. The 
algorithm for calculating non-equilibrium effects is complicated and is not repeated here. 

The thermal transpiration models were employed in stand-alone codes in order to assess im- 
pacts to system performance. Figure 4.26 shows the results of an analysis comparing the thermal 
transpiration models along the nominal reference trajectory. Here, the pressure difference between 
the nominal pressure and the pressure with the thermal transpiration effects is shown for Ports 
1 and 5. Results are shown for the Knudsen model, the abbreviated Potter model (i.e. without 
orifice effects) and the full Potter model including the orifice effects. While some difference are ap- 
parent between the Knudsen and Potter models, all show the same trend with similar magnitudes. 
These results indicate that above 850 Pa dynamic pressure, the thermal transpiration effects are 
negligible. The effects become more pronounced in the rarefied flow regime below 850 Pa dynamic 
pressure. For MEADS uncertainty analysis, the worst-case error due to transpiration is assumed 
to be 5e-4% of reading in a 3a sense. 
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Figure 4.26: Thermal Transpiration Results 
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4.1.8 System Calibration and Temperature Uncertainties 

The MEADS calibration error model is based on that described by Eq. 3.11. In summary, the 
error model consists of calibration parameter uncertainties, transducer and SEE thermocouple 
errors (noise and bias), and hysteresis and non-repeatability error estimates based on the final 
calibration. The final calibration and associated uncertainty models are described in Chapter 3. 

4.2 CFD Pressure Distribution Uncertainties 

Errors in MEADS pressure measurements due to fluid dynamic phenomena not attributable to 
MEADS hardware or software are estimated from computational and experimental results. Some 
sources of error can be estimated directly from the CFD simulations at Mars entry conditions, and 
include contributions from port location uncertainty, effects of grid refinement and topology, de- 
sign cycle changes to the OML, ablation and deformation, tile gap filler protuberances, and orifice 
induced errors. These error estimates are calculated by taking the difference between nominal and 
off-nominal solutions, and scaling to freestream conditions where applicable, or by applying func- 
tional relationships found in the literature. Each of these factors are accounted for individually and 
combined using the root sum square (RSS) method, where each MEADS port is treated seperately. 
As MEADS pressures are expressed in dimensionless form C p for scaling purposes, uncertainties in 
the pressures are expressed in normalized form A C p /C p , where the numerator is the uncertainty 
term, and the denominator is evaluated at current conditions. The RSS expression of uncertainty 
terms can be written as: 


A C p = 


^ Plocation 


+ A q 2 o„, 


^ V ablation 


+ AC Pde formation 


+ AC Lrotuberance 


+ AC? 


(4-15) 

The ability of the LAURA code to yield “correct” answers was also to be estimated by comparing 
CFD results against those obtained from wind-tunnel experimentation, and thus used to establish 
confidence in the flight CFD solutions. Significant variability in wind tunnel measurements due to 
nonuniform test section conditions were experienced. Analysis required to make the data useful for 
CFD validation purposes is beyond the scope of this work. 


4.2.1 Port Location Uncertainty 

Laser scan metrology was performed on the MSL heatshield at the Lockheed Martin Space Sys- 
tems Company in order to provide high-fidelity measurements [42] of the MEDLI/MEADS sensor 
locations and to help quantify the accuracy of heatshield-to-IMU alignment. The measurements 
indicated that the MEADS port locations are known to within approximately 1 mm, and so the 
uncertainty in their locations is quite small in comparison to other heatshield geometric features, 
e.g., discontinuities at PICA tile edges or protruding gap fillers. The effect of port location uncer- 
tainty on MEADS pressure measurement is estimated using the magnitude of pressure coefficient 
gradient extracted from flight CFD solutions, and is expressed at each port as 


A a 


| Vp| As 


Plocation 


(4.16) 


The 3cr value for port location uncertainty A C p is based on As corresponding to the metrology 
tolerance, or 1 millimeter, and is mapped over the entire nominal trajectory space. In order that 
it might be scaled to other trajectories, it is expressed in the tables normalized to the current C p , 
i.e., A Cp/ C p . 
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4.2.2 Grid Refinement and Topology, OML Design Cycle Changes 


The discretization of grids used for computation of the MSL nominal aerodatabase would be con- 
sidered rather coarse by today’s standards. After design cycle changes to the OML brought about 
as a consequence of changes to TPS requirements, a new grid of much-improved resolution was con- 
structed for the flight geometry and used to compute point checks for comparison against nominal 
aerodatabase results. The OML differences are shown in Fig. 4.27. 



Figure 4.27: Outer Mold Lines Described by Original (Grid a) and Refined (Grid b) 
CFD Grids 

Because the two grids describe different OMLs and are of differering resolution and topology, 
the effects of changing these three factors are coupled, and so are combined within one uncertainty 
term. It is important to note the limitations imposed by grid resolution and numerical stability as 
applicable to the LAURA code. On one hand, the coarse nominal aerodatabase grid is near the 
lower limit for acceptably resolved solutions; on the other hand, the updated fine grid is near the 
practical limit for computing solutions without numerical instabilities that occur with LAURA due 
to point clustering in body-tangent directions. In other words, solution accuracy will not improve 
monotonically as grid density tends to infinity. For these reasons, only the results from the original 
aerodatabase grid and the final OML fine grid are used to approximate the errors due to both grid 
effects and OML changes. The expression for the combined uncertainty term is written as 


^^POML \^P finegrid ^Pcoarsegrid I 


(4.17) 


The left hand side of Eq. (4.17 is assumed to lie at the 2c bounds, as the data were only available 
for a few flow conditions. A grid sensitivity study deemed necessary to estimate to 3c confidence 
is beyond the scope of this work. The right hand side of Eq. (4.17) is only known over a small 
portion of the trajectory space and so is implemented in a form amenable to scaling over the entire 
trajectory, i.e., A C p /C p . The maximum value of AC p /C p is used from the available data at each 
port and applied over the entirety of the trajectory space. 
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4.2.3 Ablation and Deformation 


Shape changes to the heatshield that occur due to ablation and deformation were modeled for the 
OML13F configuration. The ablation case was computed at the 07-25 trajectory peak heating 
condition, for a bounding worst-case scenario where the TPS is recessed to a maximum of approxi- 
mately 20 millimeters along the symmetry plane leeside [19]. The uncertainty term due to ablation 
is an engineering estimate defined as the absolute difference between baseline and ablated cases, 
scaled to the sum of dynamic pressure and heat load, normalized to their maximum and total 
values, respectively 

= \C Pablated — C Pbaseline \- (- b -w \ (4-18) 

^ \Qmax Wmax / 

where Q is the instantaneous heat load and Qmax is the total heat load. 

For analysis of heatshield deformation and its effect on aerodynamics, a finite element model 
subjected to 12.2 g decelerative load on the 06-05 nominal trajectory was used to modify the 
baseline grid, where predicted heatshield apex deflection is approximately 13 millimeters. As the 
current trajectory has a 3er high deceleration of 13.4 g, a scale factor is applied to the right hand 
side in order to reflect the current conditions via linear extrapolation. The expression for the 3<r 
loading term is then based upon the difference between baseline and deflected solutions, and scaled 
to the ratio of instantaneous to maximum dynamic pressure. This term is written as 


^ ^ 1 Ploading I ^Pde formed ^Pbaseline I 


Qmax 


(4.19) 


The ablation and deformation terms are only known at their corresponding peak conditions, so 
are implemented over the trajectory space using the A C v /C p form as used for the OML term. 


4.2.4 Gap Filler Protuberances 

The effect of PICA tile gap filler protuberances at maximum heating conditions on the 0801-A 
trajectory were investigated for three locations on the heatshield [43], where the nose cap case was 
deemed most relevant to MEADS performance. Heatshield flowfields with gap filler heights of 0.1 
inch, 0.2 inch, and 0.4 inch were computed at fully turbulent conditions, and curves of pressure 
along a surface streamline traversing the protrusion were extracted from the solutions, as shown in 
Fig. 4.28, where positive distance denotes downstream location. 

The radius of influence of these protrusions was shown to be on the order of 0.3 meters. For 
conditions where the MEADS port is in close proximity to a gap filler protrusion, the corresponding 
magnitude of error A C p is as large as 0.4 and dwarfs the sum of all other uncertainty terms 
considered here. In addition, the protrusion radius of influence and magnitude of error would 
increase if the flow were modeled as laminar. Uncertainty is inferred from the turbulent results and 
scaled in the same way as the ablation term, where the expression is written as 


a c. 


Pgapfille 


I ^Pgap filler ^Pbaseline 


Q 


+ 


Q 


Qmax Q 


(4.20) 


While a gap filler protuberance term would dominate the uncertainty model, the likelihood 
of its occurrence was considered low, and so the term was omitted from the model. If gap filler 
protuberances did occur, they would likely be identifiable in the pressure data as large excursions 
from the nominal model. 
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0.3 



dist (m) 


Figure 4.28: Pressure Coefficient Disturbance vs. Relative Distance from a Step Dis- 
continuity on Nose Cap. Flow is From Left to Right. 

4.2.5 Orifice Induced Errors 

The presence of a pressure orifice in a flowfield can cause errors in surface pressure measurements, 
where the errors are a function of the orifice diameter, shear stress at the wall, and the working 
fluid density and viscosity [45]. This relationship may be expressed as 

f = / ( d + ) (4.21) 

where e is the pressure error and d + is the orifice Reynold’s number, u T d/v. 

The plot of this relationship, shown in Fig. 4.29(a), shows that the left hand side becomes 
quite small (on order of 0.1) when d + is less than 100 [45]. Using the TPS stressing trajectories 
for reference conditions, the orifice Reynolds number at MEADS port locations never exceeds 
approximately 70, and so the corresponding pressure error never exceeds approximately 5 Pa. A 
plot of orifice Reynolds number is shown for the peak dynamic pressure along the most stressful 
trajectory in Fig. 4.29(b). In A C p /C p form, the peak magnitude of this error is of order l.E-4, and 
is vanishingly small for all but the lateral ports at conditions at or near peak dynamic pressure. 
For this reason, and as the curves describing this error relationship are not well established at such 
low values of orifice Reynolds number and wall shear stress, the term is omitted from the current 
model. 
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(b) d + at Peak Conditions 

Figure 4.29: Orifice Effects 
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4.2.6 Uncertainty Term Implementation 

Table 4.7 gives 3cr values for the error terms that are implemented in A C p /C p form. To obtain the 
scaled uncertainty A C p for each, they are multiplied by the base C p , which is interpolated from 
the lookup tables for the current flight condition. The terms may then be summed as in equation 
4.15. 


Table 4.7: 3a A C p Values for Uncertainty Terms Scaled Over the Trajectory Space 



Port 1 

Port 2 

Port 3 

Port 4 

Port 5 

Port 6 

Port 7 

OML + grid 

Ablation 

Deformation 

4.711e-03 

1.000e-03 

8.400e-04 

5.418e-03 

1.221e-03 

3.698e-04 

3.824e-03 

1.807e-03 

3.798e-03 

4.080e-3 

3.469e-02 

1.223e-02 

1.264e-02 

2.057e-03 

2.296e-02 

4.969e-03 

7.671e-03 

1.098e-02 

4.969e-03 

7.671e-03 

1.098e-02 


4.3 Error Analysis Results and Sensitivity Studies 

This section describes the total system performance results for both the flight and flight spare 
transducer sets. Monte Carlo analysis is utilized with the basic pressure-only MEADS algorithm 
to determine the angle of attack, sideslip, dynamic pressure, and Mach number estimation per- 
formance. Monte Carlo analysis described in this section is based on the 09-HOL-01 nominal 
trajectory. 



Figure 4.30: Pressure Measurement Errors Along Trajectory 

The total uncertainties for the MEADS flight system evaluated along the nominal trajectory 
are shown in Figure 4.30. The values of the error model parameters are listed in Table 4.8. Note 
that the error parameters are assumed to be Gaussian except where otherwise noted. 
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Table 4.8: Error Parameters 


Parameter Name 

Value 

Mensor Accuracy (a c ) (Ref. 47) 

0.001 psi 

Transducer Thermocouple Bias ) 

2.5 deg C (uniformly distributed) 

SSE Thermocouple Bias (cp, 2 ) 

1.5 deg C (uniformly distributed) 

Thermocouple Noise (a ni = a n2 ) 

0.05 deg C 

SSE Pressure Channel Noise (ay) 

10 counts 

A/D bit size 

14 

Time Stamp Error 

25 /us (uniformly distributed) 

Port Location Uncertainties 

0.04 in (uniformly distributed) 

Maximum Leak Rate 

0.047022% (uniformly distributed) 

Pressure Path Length Uncertainty (Ref. 48) 

1.41 mm 

Pressure Path Diameter Uncertainty (Ref. 48) 

0.254 mm 


The results of the Monte Carlo analysis are shown in the following figures. Figure 4.31 shows 
the total system error for the angle of attack estimate, Figure 4.32 shows the total system error 
for the angle of sideslip estimate, and Figure 4.33 shows the total system error for the dynamic 
pressure estimate. In these figures, the errors are shown only for the portion over which the MEADS 
hardware requirements pertain. Specifically the total hardware and total CFD errors are shown 
(cyan and red curves, respectively) along with the total combined error (blue). for dynamic pressures 
greater than 1750 Pa for angle of attack errors, greater than 1250 for sideslip errors, and greater 
than 850 Pa for dynamic pressure errors. These plots indicate that the MEADS hardware meets 
requirements over the entire dynamic pressure range range. 



Figure 4.31: Total Estimation Error: Angle of Attack 
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Figure 4.32: Total Estimation Error: Angle of Sideslip 



Figure 4.33: Total Estimation Error: Dynamic Pressure 

Plots of individual error sources are given in Appendix B. Here, the trajectory average errors 
for the uncertainties are tabulated in Table 4.9. These uncertainties indicate that, on average, 
the MEADS system achieves angle of attack and sideslip estimation errors less than 0.5 degrees 
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and dynamic pressure estimation errors less than 1%. The largest contributors of errors are the 
transducer performance and the OML/grid CFD errors. 


Table 4.9: Trajectory Average Errors (3cr) 

Error Model Parameter 

Angle of 

Angle of 

Dynamic 



Attack (deg) 

Sideslip (deg) 

Pressure (%) 

Transducer Errors 

Calibration Error 

0.1982 

0.1364 

0.2377 


Temperature Uncertainty 

0.0788 

0.1060 

0.1217 

SSE Errors 

Quantization 

0.1205e-ll 

0.0015e-ll 

0.7133e-ll 


Noise 

0.0488 

0.0730 

0.0971 


Temperature Uncertainty 

0.0096 

0.0082 

0.0170 

Geometry Error 

Port Location 

0.0387 

0.0301 

0.0268 


Pressure Leak 

0.0960 

0.1055 

0.1228 


Lag Uncertainty 

9.4486e-4 

1.0377e-3 

1.2138e-3 


Thermal Transpiration 

6.8808e-4 

0.7477e-3 

0.8749e-3 

Latencies 

Time Tag Error 

0.1380e-3 

0.1233e-3 

0.1560e-3 

Mechanical 

Vibration 

0.0823 

0.0966 

0.0977 

CFD 

OML/grid 

0.1966 

0.2088 

0.2376 


Ablation 

0.0426 

0.1293 

0.0400 


Deformation 

0.1294 

0.1879 

0.0941 

Total Hardware System Error 

0.2558 

0.2391 

0.3262 

Total CFD Error 


0.2392 

0.3092 

0.2586 

Total (Combined) Error 

0.3503 

0.3909 

1.3785 
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Chapter 5 

Conclusions 


This document describes the as-built MEADS error budgets for the angle of attack, angle of sideslip, 
and dynamic pressure reconstruction based on pressure measurements. The error models include 
hardware, calibration, vibration, and CFD pressure distribution uncertainties. The MEADS nom- 
inal pressure model and its corresponding CFD-based uncertainty model are presented, where 
uncertainty terms due to OML changes, grid induced error, and flight ablation and deformation 
effects are approximated using nominal aerodatabase and off-nominal CFD results. It is currently 
recommended that TPS tile gap filler protuberance effects be omitted from the uncertainty model, 
as their likelihood and magnitude are not understood. Uncertainties due to orifice induced error 
are shown to be inconsequential for MEADS at flight conditions. CFD validation and uncertainty 
models based on wind tunnel data are absent from the present work, as the variability of wind 
tunnel conditions make such models impractical to implement. Nevertheless, some improvements 
to the current uncertainty models can be suggested. These would include recomputation of the 
MSL nominal aerodatabase on the 12-GAL-01 trajectory, using the fine OML-13F grid and full 
Navier-Stokes equations. Ablation and deformation effect approximations can also be improved in 
this way. Additional uncertainty terms for code-to-code comparisons relevant to the MSL aero- 
database can also be made using different versions of LAURA, as at least two major versions have 
been released since the MSL nominal aerodatabase was computed. 

Applying Monte Carlo simulation techniques to the MEADS system model shows that the as- 
built hardware will meet Level 2 requirements over the MSL entry trajectory.The preflight modeling 
efforts described herein were validated via post-flight data processing and trajectory reconstruction 
efforts performed after MSL’s landing on Mars on August 6th, 2012. The MEDLI flight data and 
associated reconstruction activities are described in [49-52], 
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Appendix A 


Data Reduction Look-Up Tables 


This appendix contains the two-dimensional lookup tables that can be used to convert trans- 
ducer Voltage outputs into pressures (in units of psi), given the transducer and SSE temperatures. 
First, the coefficients do, ai, and a-2 are evaluated as a function of temperatures Tf (Transducer) and 
I2 (SSE), and pressure is calculated from the equation p = ao + a{V + d2V 2 where V is the Voltage 
in units on mV. Similarly, the uncertainties in each polynomial coefficient can be interpolated to 
compute a pressure uncertainty due to calibration uncertainties. 
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-128.376923 62.115385 0.155085 0.520848 -6.88512e-05 0.212737 0.0742371 0.00082291 

-128.376923 31.553846 0.149472 0.52119 -5.01051e-05 0.214292 0.0747002 0.000826397 

-128.376923 -13.030769 0.151085 0.521399 -7.63139e-05 0.213859 0.0746078 0.000826158 

-128.376923 -21.438462 0.151218 0.521272 -6.08623e-05 0.213835 0.0746218 0.000826517 
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Table A7: Transducer 026 Data Reduction Table 

Transducer SSE Coefficient Coefficient Coefficient Coefficient a o Coefficient ai Coefficient a 2 

Temperature Temperature ao ai ci 2 Variance Variance Variance 

[deg C] [deg C] [psi] [psi/rnV] [psi/rnV 2 ] [psi 2 /cr 2 ] [psi 2 /cr 2 /mV 2 ] [psi 2 /cr 2 /mV 4 ] 

-80.738462 62.115385 0.234251 0.514585 -0.000115512 0.189873 0.0680741 0.000781775 
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Appendix B 


Error Analysis Sensitivities 


Plots in this Appendix show the impact of the MEADS system errors in a single variable sense 
in order to show sensitivities to each individual error parameter. The net effects of all errors 
together are shown in Figures 4.31-4.33. Here, the MEADS hardware uncertainties are shown in 
Figures Bl-Bll and the errors due to CFD uncertainties are shown in Figures B12-B14. 




Figure Bl: Estimation Errors due to Calibration Uncertainty 
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Figure B2: Estimation Errors due to Transducer Temperature Uncertainty 
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Figure B3: Estimation Errors due to Quantization 
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Figure B4: Estimation Errors due to SSE Noise 
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Figure B5: Estimation Errors due to SSE Temperature Uncertainty) 
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Figure B6: Estimation Errors due to Port Location Uncertainty 




Figure B7: Estimation Errors due to Pressure Leak Uncertainty 
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Figure B8: Estimation Errors due to Time Tag Uncertainty 



Mach Number 


Figure B9: Estimation Errors due to Pneumatic Lag Uncertainty 
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Figure BIO: Estimation Errors due to Thermal Transpiration Uncertainty) 



Figure Bll: Estimation Errors due to Vibration 
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Figure B12: Estimation Errors due to Ablation 
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Figure B13: Estimation Errors due to Deformation 


78 



Figure B14: Estimation Errors due to OML Change 
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